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ABSTRACT 

We study the sensitivity of the processes e + e~ — > ti>i qq' at LEP2 energies on 
the non-standard trilinear gauge couplings (TGC), using the optimal observables 
method. All relevant leading logarithmic corrections to the tree-order cross section, 
as well as experimental resolution effects have been studied. Taking into account 
correlations among the different TGC parameters we show that the limits on the 
TGC can reach the level of 0.15 (lsd) at 161 GeV with 100 pb _1 , a challenge for the 
first LEP2 phase. At higher energies this can be improved drastically, reaching the 
level of 0.02 (lsd). 
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One of the two most important measurements at LEP2 energies will be the determination of 
the trilinear vector boson couplings jl], ^, 0], a characteristic manifestation of the underlying 
non-Abelian symmetry of elementary particle interactions ||. 

In order to study the trilinear boson couplings we need a parametrization of these interactions 
that goes beyond the Standard Model . There are, of course, several possibilities to accomplish 
this task, but we are going to restrict ourselves to the most economical one by considering C— 
and P— conserving interactions. The relevant interaction Lagrangian is usually written in the 
following form ||: 

Ctgc = E e9v(V^W-^W+ -V^W+^W- + k v V imu W + »W-») 

V=i,Z 

+ e^A, p W + ^W-» + ectg6 w ^Z, p W + ^W-» (1) 

where 

v^ = d^ v -d v v p , w% = d^wt - d v w± . 

is the VF-boson field, and = 1, gz = ctg 9 W , k 7 = kz = 1 and A 7 = Xz = at tree 
order in the Standard Model . It is more convenient to express the different couplings in terms 
of their deviations from the Standard Model values. For this we define the following deviation 
parameters^]: 

= gz - ctg 8 W x 1 = k 1 -1 x z = (k z - l)(ctg 9 W + 5 Z ) ■ (2) 

It is worth while to note that the interaction Lagrangian becomes linear with respect to the 
above parameters (including also A 7 and Xz)- 

During the last years, considerable progress has been achieved concerning the understanding 
of the physics underlying the non-standard boson self-couplings. As Gounaris and Renard || 
showed, the deviations from the Standard Model couplings can be parametrized in a manifestly 
gauge-invariant (but still non-renormalizable) way, by considering gauge- invariant operators 
involving higher-dimensional interactions among gauge bosons and Higgs field. These operators 
will be scaled by an unknown parameter A^p, which might be understood as the characteristic 
scale of New Physics effects. In order to describe all five C— and P— conserving couplings 
introduced in Eq.(|l]) we need operators with dimension up to eight. On the other hand restricting 
ourselves to SU(2)l x [/(l)y-invariant operators with dimension up to six, which are the lowest 
order ones in 1/Ajvp expansion, we can have the following list Q: 

Ob* = B^{D^)\D v <b) 
Ow* = (D^ t-W^(D u $) 

Ow = ^(W^xWV-W^ (3) 
where Tj = \(Ji (cij are the Pauli matrices), 

B^u = d p B v — d v B p 

where B^ is the U(l)y gauge field, 

W.u = duW v -d u W„- gW, x W v 
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Figure 1: Single-resonant graph where TGC are contributing, 
where W are the SU(2)l gauge fields and 

$= ( ^(v + H + iitP) ) 
is the Higgs doublet. The covariant derivative is given, as usual, by 

Dp = df. + i gr ■ - i g'B^ 

and e = g sin W = g' cos 6 W . 

The interaction Lagrangian can be written now as 

Ltgc = ig -^-Vb$ + ig — — <-V* + g—^-Uw (4) 

m W m w m w 

where the relations between aw$, ob$, aw and the deviation parameters of Eq.(||) are given 
by: 

$z = a-w$/ (s'm9 w cos 9 W ) x 7 = ag$ + aw<s> A 7 = aw 

xz = — tan 9 w x~/ Xz = A 7 (5) 

In order to study the effect of TGC, one traditionally considered the reaction e + e~ — ► 
VF + VF~, taking into account the subsequent decay of the W's in a four-fermion final state. 
These final states can be classified in three categories, namely the 'leptonic' l\VgJG±y^ the 
'semileptonic' tvi qq' and the hadronic qiq'i&q^ (q an d q' refer to up- and down-type quarks 
respectively). Semileptonic seems to be the most favoured channel || for studying TGC, since 
it contains the maximum kinematical information, taking into account that charge-flavour iden- 
tification in four jet channel is rather inefficient and the cross section for the leptonic mode 
is suppressed. It is the goal of the present paper to study the effect of TGC, in their three 
parameter version Eq.(|j), in the processes 

e + e~ — ► £Pi qq* (6) 

where I is an electron or a muon, at LEP2 energies, based on the four-fermion Monte Carlo 
generator ERATO [10|, O]. The final state rv T qq' will not be considered here due to the 
special difficulties to identify r's in this environment. 

From the point of view of the four-fermion production, e + e~ — ► W + W~ calculations are 
equivalent to the narrow width approximation, Tw — * 0, which requires that both Ws are 
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kinematically allowed to be on-shell. This last requirement implies of course that the overall 
energy should be above the 2mw threshold and therefore it is not adequate for the first phase 
of LEP2, where yfs = 161 GeV. Moreover in e + e~ — ► the invariant masses of the 

produced fermion pairs, tvg and qcf , are identical to the W mass, which requires experimental 
selection algorithms, whose efficiency is some times questionable, and more importantly their 
interplay with TGC studies can not be estimated, unless a full four-fermion calculation is used. 
On the other hand there are nowadays widely available four-fermion calculations (9|, 12, 10 1 where 



TGC are included beyond the narrow width approximation. This enable us to study TGC not 
only in the double-resonant (CC3) graphs contributing to e + e~ — > W + W~ but also in the 
single-resonant ones shown in Fig.[I], which moreover become dominant at very high energies. 

In the actual calculations presented in this paper we have used the ERATO Monte-Carlo 
generator. A detailed description of ERATO can be found in references || [l0|, |llj. The basic 
ingredients of the calculation are the following: 

1. Exact tree-order matrix elements for the processes e + e~ — > tvg qq' , based on helicity 
amplitude calculations 0, § , including all TGC interactions described by Eq. (||) . 

2. Phase-space generation algorithm based on a multi-channel Monte Carlo approach, includ- 
ing weight optimization [13]. 



3. Treatment of unstable particle (W and Z) width consistent with gauge-invariance and 
high-energy unitarity |l4|, [H| . 



4. Initial state radiation (ISR) based on the structure function approach] 16, [T(J, including 
soft-photon exponentiation as well as leading logarithmic (LL) contributions up to order 
0{a 2 ). 

5. Coulomb correction^ to the double resonant (CC3) graphs [16]. 

Moreover in order to avoid matrix element singularities as well as to be as close as possible 
to the experimental situation, the 'canonical cuts' of reference [10] have been applied: 



175° > 9 t , Ojet > 5°, Ei > 1 GeV , E jet > 3 GeV and m qA , > 5 GeV (7) 

The input Standard Model parameters we have used are as follows: 

M w = 80.23 GeV, T w = 2.033 GeV, M z = 91.188 GeV, T z = 2.4974 GeV, 

sm 2 9 w = 0.23103 and a = 1/128.07 (8) 

As far as the ISR structure function is concerned the value a = 1/137.036 is used. 

Most of the techniques used up to now to determine the sensitivity of a given reaction on 
the TGC have been based on a Maximum-Likelihood fit to Monte-Carlo generated four-fermion 



data. In its most advanced version 17], which is capable of dealing with low statistics data 



without any binning procedure, one maximizes the log-likelihood function: 

ln£ M L = ^2ln p(Qi, a) , (9) 



1 For a detailed description see ref.Jlf]], pp 117-119. 
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where the sum is running over the event sample, fij represents the collection of the independent 
kinematical variables describing the i-th event, a = (aw$, aB$, aw) is the vector whose compo- 
nents are the three TGC parameters, p(ftj, a) is the probability to find an event at such a phase 
space point, and is defined by: 



p(ft,a) 



1 da 
a dfl 



(10) 



with 



a 



da 



dQ . 



Although this is the most direct method to analyse experimental data, one has to calculate 
the full differential cross section for a large number of TGC parameter values, which is rather 
time-consuming. 

On the other hand, since the interaction Lagrangian is linear with respect to TGC, one can 
easily write the differential cross section in the form 



da 

dtt 



co 



(ft) 



(11) 



and of course in a similar way the total cross section is written as 



(12) 



where hatted c's are integrals of unhatted ones over the phase space. It is now straightforward 
to verify that the probability has an expansion 



Cl4P) _ Cl£ 

co (ft) c 



+ higher orders 



where 



1 da 

Po{n) = adn 



(13) 



(14) 



a=0 



is the Standard Model distribution. 

Assuming now that the data are well described by the Standard Model (a = 0) distribution 
it is easy to calculate the so-called information matrix 18 j, which determines the sensitivity of 
the fit on the parameters a, and is defined by: 



E 



-S- In C M l 

da; 



_d_ 

da-i 



ln£ 



ML 



(15) 



-E 



d d 
den dcij 



ln£ 



ML 



where 



N 

E[A] = J n{rfftiPo(fti)M(fti>->ftiv) 

i=l 



(16) 



represents the mean value of a function A. To the lowest order one has that J« = NBij where 



Bi 



(^^)o-(^)o(^)o 
CO c c c 



(17) 
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and 

(A) = I dn Po (n)A(n) . (is) 



Using now Eq.(|l3|) one has that the set of the phase-space functions: 



a, = ■ (M) 

Cq{U) 

are the optimal observables []l9| , whose mean values and covariance matrix will determine the 
sensitivity on the TGC parameters, and it is equivalent, to the lowest order, to the log-likelihood 
method. More specifically we have that 

{O i ) = (O i ) + J £a j B ij . 

j 

The estimator is now given by 

j 

and the corresponding covariance matrix 

V(a) = 1 B- 1 ■ V(0) ■ B- 1 
where V{0) is the covariance matrix of the Oi defined by 

Vip)i3 = {p i O j )-{O i )(O j ) 

which in the ideal case, where measured distributions are identical to the Standard Model 
predictions, is given by V(0) = B. 

One-dimensional investigations assuming that all but one of the a's are non- vanishing, result 
in parameter errors (lsd) given by 

Sa - - vm ■ m 

On the other hand taking into account the correlations one has to diagonalize the covariance 
matrix V(a) and therefore 

Sa -=7m (21> 

where Aj are the eigenvalues of the matrix B and the parameters af are defined by 

af = ei-a . (22) 

where e*j are the corresponding eigenvectors of B. 

In all calculations presented in this paper, ./V is taken to be the predicted number of events 
in the Standard Model , defined as 

N = 4 L a 

where a is the corresponding total cross section, L is the integrated luminosity and the factor 
4 takes into account the four equivalent channels (all described by the same matrix elements), 
i.e. e + e~ — > Ivi qq' stands for e + e~ — tv^ ud, £Pi cs, tvg du and £*V£ sc. 
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At this point it is worth to emphasize that expanding the log-likelihood function to higher 
orders in eij will allow us to identify observables which probe the behaviour of log Cm L up 
to the corresponding order. Therefore there is no limitation, in principle, to establish a one- 
to-one mapping between Maximum-Likelihood and optimal observables. On the other hand 
one can argue that since we have truncated our interaction lagrangian to the order l/h? NP , 
consistency would require that each observable should also be trusted up to this order and 
therefore only linear terms should be kept. In any case optimal observable and Maximum- 
Likelihood become identical in the high statistics limit, N — > oo, and therefore we are well 
justified to use optimal observables to estimate the errors on TGC parameters. Moreover, the 
optimal observables method can easily be extended |l9|], when the variation of the total event 
rate is taken into account, corresponding to the so-called Extended Maximum Likelihood (EML) 
method. It should be noted, however, that in the case where some of the en terms vanish, or 
<(|) 2 > " & ~ 0(1/ y/N) x ((aa) _ leading and next-to-leading terms become 

equally important and an analysis based on the optimal observables of Eq.((l^) breaks down. 
Fortunately, such a 'pathological' case is not observed at LEP2 energies. 

Finally, from the point of view of a weighted Monte-Carlo approach, which is frequently used 
in the phenomenological analyses, optimal observables offer a very efficient, fast and economic 
way to estimate the sensitivity of a given process on the TGC (or any kind of 'deviation' 
parameters). The only thing one has to do is to calculate the correlation matrix Bij. This 
should be contrasted with log-likelihood methods where either one has to repeat the calculation 
for a large set of TGC parameter values which adequately covers the parameter space, or to 
store information on the coefficients appearing in the Eq,(|ll~|) for each generated 'event'. Of 
course from the point of view of 'unweighted' or real data analysis these arguments do not apply 
any more, since now local (event-by-event) information is any way available. In this case log- 
likelihood would be preferable over optimal observables, since it is not restricted to the linear 
terms. 





161 


176 


190 


205 


e 


0.456 0.180 0.087 
0.180 0.173 0.008 
0.087 0.008 0.117 


0.587 0.086 0.238 
0.086 0.087 0.014 
0.238 0.014 0.311 


0.973 0.115 0.420 
0.115 0.135 0.024 
0.420 0.024 0.605 


1.400 0.163 0.601 
0.163 0.209 0.036 
0.601 0.036 0.969 


124.6(6) 


489(2) 


591(2) 


621(2) 


A 4 


0.271 0.013 0.092 
0.013 0.019 0.004 
0.092 0.004 0.093 


0.547 0.033 0.253 
0.033 0.037 0.014 
0.253 0.014 0.313 


0.956 0.063 0.447 
0.063 0.084 0.028 
0.447 0.028 0.620 


1.435 0.111 0.660 
0.111 0.151 0.044 
0.660 0.044 1.036 


123.1(5) 


476(2) 


567(2) 


585(2) 



Table 1: The correlation matrix for e and \x channels at LEP2 energies for TGC parameters 
a = (aw®, a_B$) clw)- Also shown are the cross sections as well as their Monte Carlo errors in 
femtobarns. 
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In table [l] we present the results for the correlation matrix, Bij , for centre of mass energies 
161, 176, 190 and 205 GeV. The total cross sections are also presented. As is evident from 
this table, the correlation between the different a's are rather strong, which suggests that an 
analysis taking into account these correlation is indispensable. As is easily seen from this table, 
the correlations are much stronger between aw$ and aw as compared with clb$, which in general 
gives the smallest contribution. 

An other very interesting result, is that electron and muon channels exhibit different be- 
haviour depending on the energy. The difference is more prominent at the threshold region, 
y / i=161Gev. This is due to the contribution from single-resonant VT-production graph Fig.[l], 
which becomes relatively more important at threshold. On the other hand at higher energies, 
the dominance of the double-resonant graphs (CC3) restores universality among electrons and 
muons, and the correlation matrix is almost identical for both channels. It is worth to empha- 
size that differences in the correlation matrix (~ 100%) among e and /u channels are much more 
important than the differences in the total cross section (~ 1%), reflecting the fact that the 
graph of Fig.|l] contributes mainly to the shape of the differential cross section rather than to its 
overall normalization. 





e 




161 


0.559 (0.889, 0.419, 0.183) 
0.125 (0.128,-0.613, 0.779) 
0.062 (0.438,-0.670,-0.599) 


0.311 (0.920, 0.046,-0.390) 
0.054 (0.388, 0.039, 0.921) 
0.018 (0.058,-0.998,-0.018) 


176 


0.734 (0.863, 0.125, 0.489) 
0.183 (0.443, 0.274,-0.853) 
0.068 (0.241,-0.953,-0.181) 


0.711 (0.841, 0.052, 0.537) 
0.152 (0.535, 0.051,-0.843) 
0.035 (0.072,-0.997,-0.015) 


190 


1.258 (0.835, 0.097, 0.540) 
0.339 (0.517, 0.193,-0.833) 
0.116 (0.186,-0.976,-0.111) 


1.269 (0.821, 0.057, 0.568) 
0.311 (0.565, 0.055,-0.823) 
0.080 (0.078,-0.997,-0.013) 


205 


1.838 (0.817, 0.094, 0.569) 
0.557 (0.549, 0.173,-0.817) 
0.183 (0.175,-0.980,-0.089) 


1.932 (0.802, 0.065, 0.594) 
0.548 (0.589, 0.075,-0.804) 
0.141 (0.097,-0.995,-0.023) 



Table 2: The eigenvalues and the corresponding eigenvectors of the correlation matrices given 
in table H. 

In table § we show the eigenvalues of the correlation matrix, as well as the corresponding 
eigenvectors. These eigenvectors define directions in the three-parameter space, which are uncor- 
rected and therefore parameter errors can be safely extracted. As we can see the two dominant 
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eigenvalues correspond to directions in the three-parameter space related predominantly to a\y$ 
and aw, whereas the lowest eigenvalue is always related to 05$. This ob<j> suppression is due to 
the fact that in the double resonant (CC3) graphs relative cancellations between terms propor- 
tional to x-f and xz take place^j. An exception to this pattern is due to the electron channel at 
161 GeV, where the mixing is much more important, reflecting the contribution from the graph 
of Fig.|l|. This is easily understood, since in the amplitude corresponding to this graph the rela- 
tive cancellations between x 7 and xz are now destroyed due to the the presence of the i-channel 
photon propagator, which means that the f-channel single-resonant graph receives contributions 
only from x 7 and A 7 whereas contributions from xz, Sz and Xz are suppressed. 

Finally in table ^ one-standard-deviation errors are presented by combining e and /x in the 
following wayQ 



J v "U a (e) + U (a 1 ) ^ « CT (e) + <r(/0 ' 
and 

N = 4L(a (e) . 

In all our studies, L is taken to be 100 pb" 1 at 161 GeV and 500 pb _1 at all higher LEP2 ener- 
gies. In ld-case, corresponding to 'one-dimensional log-likelihood fit', errors on a\y§, aw, clb^ 
are shown. For 3d-cases, the errors correspond to the linear combinations defined in Eq.(22). 
Nevertheless in table || we kept the notation of aw$, aw and 05$ for 3d-cases, since the cor- 
responding eigenvector, a w § for instance, is mainly composed by ajy$ and so on for the other 
TGC parameters. It is rather evident that correlations, though important, do not alter dramat- 
ically the estimate of the errors based on simple 'one-dimensional' investigations, and therefore 
Id limits are still useful. 





aw§ Q>w a B§ 
161 


aw& aw 
176 


aw® aw o_b$ 
190 


aw§ o>w a B& 
205 


Id 


0.16 0.31 0.32 


0.030 0.041 0.091 


0.021 0.026 0.062 


0.017 0.020 0.048 


3d 


0.15 0.32 0.44 


0.027 0.056 0.098 


0.018 0.036 0.066 


0.015 0.027 0.050 


3d' 


0.18 0.37 0.50 


0.031 0.064 0.110 


0.021 0.041 0.073 


0.017 0.030 0.055 



Table 3: One standard deviation errors on TGC parameters. Id means that the correlation 
matrix is assumed diagonal where for 3d and 3d' all correlations have been taken into account. 
Finally in 3d' 'resolution' effects are also included as described in the text. 



In real life, detector resolution effects and reconstruction algorithm, tend to distort the 
'observed' distributions as compared with the 'theoretical' predictions. This can in principle, 
mimic the effects of TGC, and therefore should be studied on a generator level. This study 
requires a good knowledge of the specific detector and should be carried out in detail by our 
expiremental colleagues. Nevertheless, we can estimate the order of magnitude of these effects 

2 that is cancellations between 7 and Z-exchange, s-channel graphs. 

3 which becomes exact when (ci,i/co)(cij/co) <C {(ci,j/co) (ci,j/co)), as is indeed the case in our calculations. 
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based on the following very general assumption, that the distorted-'observed' distribution p e ff 
is expressed as a convolution integral given by: 

Pef fi^meas > fl) — J PiS^truet Si^ltrueT^meas) ddf rue (^3) 

where p(0, true , a) is the theoretically predicted distribution and S(p, true , Q mea s) is a 'resolution' 
function, giving the probability that an event 'measured' at phase space point £l m eas is coming 
from a 'true' event at a phase space point Qtrue- 

In order to simulate in a simple way the 'resolution' function we have implemented the 
following algorithm: 

• To each generated event we calculate the nine kinematical variables, in the lab frame, 
Ei,8i,(j)i, i = 1,2,3, which correspond to the 'visible' particles, namely the lepton £ and 
the two jets q, q' . 

• To each of the above variables we assign a new one, 

9', 4>^4> and E -> E' 

following the rule that primed variables are normally distributed around the original ones^ 
with variances given by: 

<T Et = 0-03 Et, a 6e = 0.0002, a 4>l = 0.0002, 
a Ejei = 0.15 E jeU a 6jet = 0.03, a^ jet = 0.04 . 

• Then, this artificially generated phase space point is subjected to a 'kinematical constraint 
fit' algorithm, by including a 'missing energy' (neutrino) four-momentum in such a away 
that 

4 4 

^2 Ei = 2 E beam and p f = 

i=l i=l 

This is achieved by a local (event by event) rescaling of the energies of the two jets, 
which are expected to be the worst measured quantities. This rescaling is defined by the 
minimization of the following x 2 function: 

X ~ 2^ a 2 

jets 1 

Note that, although the original phase space point is generated taking into account initial 
state radiation, in the 'reconstructed' one we have artificially neglected it and the total 
energy is normalized to the overall constant beam energy. 

• Finally the usual cuts of Eq.(0) are applied, including now a special cut on the missing 
far vector, p^r > 15 GeV. This supplementary cut is applied in order to minimize the 
background coming from the process e + e~ — > e + e~ qq, where the e + is assumed to be 

4 with the exception of jet energies, where (E') = 0.85 E, in order to take into account energy losses in jet 
reconstruction. 
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lost into the beam, which leads to the same final state e~ + 2jets. Using ERATO we 
have checked that the reconstructed for in the background process is well below 15 GeV, 
whereas the ]/>t distribution of our 'signal' processes e + e~ — ► tvi qq has a negligible tail 
below 15 GeV. 

Using now the kinematical information from the above algorithm we can calculate the new 
correlation matrix, B^, using Eq,(|23]), and from it the new information matrix given by 

jnew j _ j— 1 _ j 

The lsd-erros on TGC parameters calculated through I new are given in table |3|, under the 
3d'-case. Comparing with the 3d case, we see that the effect is of the order of 10-20% on the 
parameter errors, which is at an acceptable level. We have checked that reasonable input values 
for the variances cte, &e an d cr^ give more or less the same results. The same is still true if 
we replace the polar representation of the observed momenta by the cartesian one. This means 
that resolution effects are unlikely to destroy the measurement of TGC. Moreover it is worth to 
emphasize that 'unfolding' techniques can also be used in order to minimize this effect, relying 
on a good knowledge of the detector as well as of the physics included in the event generator 
programme. 

We conclude by summarizing the results of our study: 

1. Sensitivities on TGC at LEP2 are of the level of 0.15-0.5 (lsd) at 161 GeV, making the 
threshold phase of LEP2 the best up to date world measurement of TGC. Higher energies 
will drastically improve these limits, almost by an order of magnitude, reaching the level 
of 0.02. Taking into account information from other channels, such as the four-jet one, 
further improvement is expected. 

2. The electron channel due to single-resonant W contributions is rather important at 161 
GeV and special studies are needed on the experimental side, in order to maximize se- 
lection efficiencies in this channel. Finite-width, ISR and Coulomb corrections are also 
indispensable at 161 GeV. 

3. Resolution effects are unlikely to drastically alter this picture, leading to a mild loosening 
of the sensitivity limits by a 10-20% factor. 

As a postscript of our study we present in Fig.|| the distributions of the optimal observable 
O awi , at 161GeV. The solid line corresponds to the ideal case, where the measured distribution 
is assumed to be identical with the predicted one, whereas the dashed line is calculated by the 



convolution integral Eq.(23), based on our naive 'resolution' model. Also shown are the two- 
dimensional scatter plots, between the generated and 'measured' values of O aw<s> . Although our 
'resolution' model is very simplistic, it suggests that optimal observables can be used in the 
experimental analysis of four-fermion data. Optimal observables are the 'phase-space' variables 
exhibiting the maximum sensitivity on TGC and therefore they effectively reduce the number 
of independent phase-space variables from eight to three, which is rather appreciable from the 
experimental point of view. On the other hand their sensitivity on 'resolution' effects seems to 
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be rather mild. Nevertheless their usefulness has to be verified by a detailed analysis at the 
experimental level. 
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Figure 2: The distribution of the optimal observable corresponding to the parameter aw<s> at 
161 GeV, for electron (a) and muon (c) channels. Also shown are the scatter plots (b) and (d). 
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